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Abstract. Despite the significance of Type la supernovae (SNela) in many fields in 
astrophysics, SNela lack a theoretical explanation. The standard scenarios involve ther- 
monuclear explosions of carbon/oxygen white dwarfs approaching the Chandrasekhar 
mass; either by accretion from a companion or by a merger of two white dwarfs. We 
investigate the contribution from both channels to the SNIa rate with the binary popula- 
tion synthesis (BPS) code SeBa in order to constrain binary processes such as the mass 
retention efficiency of WD accretion and common envelope evolution. We determine 
the theoretical rates and delay time distribution of SNIa progenitors and in particular 
study how assumptions affect the predicted rates. 



1. Method 



Binary population synthesis (BPS) codes are very useful tools to study the evolution 
of binary stars and the processes that govern them. They enable the simulation of 
large numbers of binary systems and with that the study of populations of e.g. no- 
vae, X-ray binaries and Type la supernova (SNela). BPS codes are particularly useful 
to study trends in stellar populations and the effects that various important (and un- 
certain) phases have on the binary evolution. Here we investigate the population of 
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SNIa progenitors from two possib l e evolutionary channe l s using the BPS cod e SeBa3 
dPortegies Zwart & Verbuntl 11996); iNelemans et al1l200ll; iToonen et all 1201 2h . In the 
single degenerate (SD) channel (IWhelan & ibenlll973l) a carbon/oxygen white dwarf 
(WD) accretes from a hydrogen-rich companion to approach the Chandrasekhar mass. 
The do uble degenerate (DP) channel des cribes the merger of two carbon/oxygen white 
dwarfs (|Iben & Tutukovlll984l : IWebbinklll984l) . Other channels are a.o. accreting WDs 
from helium-rich non-degenerate companions ( Wang et all I2009L see e.g.) and the 
merg er between a WD and the core of an asymptotic giant branch star dKashi & Sokei 
2011). We study the SNIa delay time distribution (DTD), where the delay time is the 
time between the formation of the binary system and the SNIa event. In a simulation of 
a single burst of star formation the DTD gives the SNIa rate as a function of time after 
the starburst. The DTD is linked to the nuclear timescales of the progenitors and the 
binary evolution timescales up to the merger. 



2. Double degenerate channel 

We set out to predict SNIa rates for the DD channel, with the additional constraint that 
our model corresponds well to the observed population of close double WDs (CDWDs) 
- of all masses and compositions. Even though there are no certain DD SNIa progenitors 
among the observed CDWDs e.g iMarshl (1201 lh . the DD SNela and observed CDWD 



progenitors have gone through similar evolutionary paths and are affected by the same 
binary and stellar processes. 

The CDWDs are believed to engage in at least one phase of common envelope 
(CE) evolution. In spite of the importance of the CE phase, it remains poorly under- 
stood. Several prescriptions for CE evolution have been proposed. The a formalism 



(Webbink 1984j) is ba sed on the conservation of orbital energy and the y-formalism 
(INelemans et al.ll2000l) is based on the conservation of angular momentum. In model 



aa the a-formalism is used to determine the outcome of every CE. For model ya the 
y-prescription is applied unless the binar y contains a c ompact object or the CE is trig- 
gered by the Darwin-Riemann instability (lDarwinlll879h . The ya-model reproduces the 
mass ratio distribution of the observed DWDs best, see Fig. [T| 

The DTD of model ya and model aa are similar in showing strong declines with 
time, approximately oc I ft, see Fig. [2^. This is expected when the delay time is dom- 
inated by the merger time of a double WDs as a result of gravitational wave emission. 
Fig. [2k shows that these mergers are expected to take place in young as well as old pop- 
ulations. In absolute numbers the DTD of model ya is similar to model aa as well; the 
time-integrated rates are 2.0-10 _4 M- 1 for model ya resp. 3.3- lO" 4 ^ 1 fo r model aa. 

The recovered DTD from a compilation of observations and method by Maoz & Mannucci 



(201 1) has a similar shape as our synthetic DTDs, bu t has a normalizati on that i s a fac - 



tor -10 higher. The integrated rate from iMaoz et al.l d2010h : iMaoz & Mannucci d201ll) 



is 18 - 29 • 10 -4 Mq 1 resp. > 34 • 10~ 4 Mq 1 based on galaxy cluster measurements and 
cluster iron abundances respectively. However, r ecent measurements i n field galaxies 
have shown lower rates; 0.6 - 1.0 • 10" 4 M^ 1 by Perr ett & et all (l2012h as rescaled by 



1 SeBa is incorporated in the Astrophysics Multipurpose Software Environment, or AMUSE. This is a com- 
ponent library with a homoge neous interface structure, and can be downloaded for free at amusecode . org 
iPortegies Zwart & et al.ll2009l) . 
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Figure 1 . Simulated distribution of population of visible double white dwarfs as 
a function of orbital period and mass ratio, where mass ratio is defined as the mass 
of the brighter white dwarf divided by that of the dimmer white dwarf. Left model 
ya is used, on the right model aa. The intensity of the grey scale corresponds to 
the density of objects on a line ar scale. The same grey scale is used for both plots. 
Observed binary white dwarfs (Marsh 2 011b are overplotted with filled circles. 



iMaoz et all (EoT 
by iGraur & Mao 
lower limit) 



13 + 1.5 • 1(T 4 M~ x bvlMaoz et all d2012l) and 5 ± 2 • 10" 4 M' 



(120121) (however IGraur & Maozl (12012) state this value is likely an 
At this moment it is unclear if the different observed integrated rates are 
due to s ystematic effects or if there is a real enhancement of SNela in cluster galaxies. 
See also Mao z et al.1 ((2012) for a discussion. 



The normalisation of the synthetic SNIa rates is dependent on assumptions in the 
simulation; for example the stellar winds, CE evolution, assumed binary fraction and 
the initial distribution of masses and orbital parameters. However, preliminary results 
show that the integrated rates are not affected by factors as large as ~ 10 (see also Claeys 



et al. i n prep.). Therefore the synthetic rates are inconsistent with the rate of IMaoz et al. 



(201 1) and the main contribution to the SNIa rate has to come from other channels. On 
the other hand, the synthetic rates are consistent with the lower limit of the SNIa rate in 
field galaxies. Future observations will help to solve this enigma. 



3. Single degenerate channel 

The single degenerate scenario explains in a very natural way the uniformity in the 
SNIa lightcurves, as all SD events occur when the WD approaches the Chandrasekhar 
mass. The theoretical SD rates predicted by different BPS codes, however, vary over 
four orders of magnitu de and do not agree with each other or with observational data 
(Nel emans et al.ll2012h . The exact origin of these differences remains unclear so far. 

A crucial phase in the evolution of SD systems is the phase of mass transfer and 
accretion onto the WD. The efficiency of mass accretion, hereafter retention efficiency 
r], is poorly understood because of processes such as novae and stable burning^. If 
rj = no mass is retained by the WD, when n = 1 the WD accretes matter conser- 



2 Note that the retention efficiency is of lit tle importance to the D P channel and that the predicted DTDs 
from different BPS codes are comparable (Nelem ans et al.l2012l) . 
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Figure 2. SNIa rate per yr per M formed stellar mass as a function of delay time, 
a) On the left rates from the DD channel. The delay time is in Gyr. The solid line 
represents a simulation using model ya for the common envelope evolution, dashed 
aa. Overplotted with bl ack circles are the observed values of the SNIa rate (see 
M aoz & Mannuccill201 ll for a review), b) On the right rates from the SD channel. 
The delay time is in yr (logarithmic). Note that no SNela where found when the 
retention efficiency of Yungelson (2010) are used. 



vatively. In BPS codes different prescriptions dNomoto et al.1120071; iRuiter et ail 12009 



Yungelsonll20 1 0|) are used for the retention efficiency. The retention effici ency of hydro- 



gen ac cretion that is burned into helium is shown in Fig. |3^, soli d 77hi for Nomot o et al. 
(120071) and dashed r/m as assumed by IRuiter et all d2009l) and lYungelsonl d2010h . At 
mass transfer rates M < 10 _7 M©yr _1 novae occur taking away a large fraction of the 
accreted mass. At mass transfer rates of a few times 10~ 7 M o yr _1 hydrogen burns sta- 
bly into helium on the surface of the WD. At higher mass transfer rates, strong winds 
develop from the surface of the WD. In th e wind regime rim > ?7H2 because a stripping 
effect of the donor by the wind of the WD ( Hachisu et al.ll2008 ) is taken into account in 
?7hi- After a fraction of the hydrogen-rich matter is burned into helium-rich material, a 
fraction of the helium-rich material is burned into carbon-rich matter. For the efficiency 
of helium retention ?7He two pres criptions exist as we ll. T he dot-dashed line in Fig. [3^ 
shows 77He = ?7Hei as assumed by iNomoto et al.1 d2007l) and lRuiter et al.1 d2009h and 77^2 
assumed by Yun gelsonl d2010l) . 

The total retention efficiencies 77 tot = ?7h(^h) • ?7He(^He) after hydrogen and he- 
lium burning is shown in Fig. [3]). This figure shows a large v ariety in the assumed 
retention efficienci es, where the mo st optimistic case is given by iNomoto et al.l d2007l) 
and pessimistic by Yun gelsonl (2010). We have used the BPS code SeBa to simulate the 
SNIa rates in the SD channel for each retention efficiency. The simulated SNIa rates are 
significantly affected by which prescription for the mass retention efficiency is used. In 
Fig. |2j) the resulting DTDs are shown assuming a y-CE. The integrated rates vary be- 
tween 5.8- 10" 5 M" 1 , 1 9- 10" 5 M" 1 and an u p per limit of 1 • 10~ 7 M" 1 when the re t ention 
efficiency is used as in lNomoto et all d2007h . Ruiter et all d2009h and lYungelsonl d2010h 
respectively. This study shows that depending on the assumptions for the retention ef- 
ficiency the co ntribution to the S NIa rate from the SD ch annel can be negl igible (for 
lYungelsonlboiOl) to moderate (for lNomoto et al.ll2007l and IRuiter et al.ll2009l where the 
difference in the integrated rates is of the order of a few). Predictions from BPS codes 
on the contribution from the SD channel to the SNIa rate are dependent on many as- 
sumptions e.g. CE evolution, but is crucially dependent on the retention efficiency. For 
more information about this study, see Bours, Toonen, Nelemans in prep. 
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Figure 3. The mass retention efficiency of WD hydrogen accretion as a func- 
tion of the mass transfer rate. rj lot represents the fraction of the transferred mass 
that is retained by the WD. In the figure we have assumed the whi te dwarf to be 
of one solar mass 



a) Left the hydrogen retention effic iency h hi (Nomoto et al, 
2007|); lHachisu et ail d2008l) - solid line) & rj m dRuiter et all d2009t). lYungelson 



2010) - dashed l i ne) an d the helium retention effici ency ?7He i (Hachisu et al. 



iKato & Hachisul i 19991) - dot-dashed line) & 77^2 (llben & Tutukovl (119961) - dotte d 
line) are shown, b) Right the total reten tion efficiency based on lNomoto et al.1 {2007); 
iRuiter et"aD (120091) : lYungelsor] (1201 0) are shown resulting from different combina- 
tions of the hydrogen and helium retention efficiencies. 



4. A comparison of binary population synthesis codes 



In the previous sections we have discussed the effect on the SNIa rate of two phases in 
binary evolution that are not well understood; the CE and mass accretion onto WDs. 
However, we are not able to reproduce the same DTDs and progenitors as other BPS 
codes when assuming their assumptions for the CE and rj. To understand the differences 
in the predictions of the various BPS codes, we started a collaboration to compare four 
BPS codes. The goal is to investigate whether differences in the simulated populations 
are due to numerical effects, or whether they can be explained by differences in the 
input physics. Regarding the latter, we aim to identify and qualify which assumptions 

in BPS codes affect the results most, and hence should be studied in more detail. 

The codes involved are the Binary _c cod e dlzzard & Toutll2004l; llzzard et al.|[ 2006. 



2009; IClaevs et al.l201 1 ), the Brussels code (j Pe Donder & Vanbeveren 2004; Mennekens et al 



2010). SeBa (iPortegies Zwart & Verbuntlll996l;lNelemans et al.l2 00ll;lToone n et al.l2012l) 
and Startrack (Belczynski et al.l l200l l2008TlRuiter et al.ll2009l) . The comparison fo- 
cuses on the evolution of low- and intermediate-mass binaries containing one or more 
WDs. For this project the four codes have equalized their assumptions as much as pos- 
sible, even when these assumptions are not assumed to be realistic. These assumptions 
are: 

• The same initial distribution of primary mass, mass ratio and separations; 

• The same CE prescription and efficiency; 

• Stable mass transfer is conservative and to all types of stars; 

• No tidal effects, magnetic braking or wind accretion. 

Preliminary conclusions are that when the input assumptions are equalized, the differ- 
ent BPS codes give similar populations. There are small differences, however they are 
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not caused by numerical differences, but can be explained by inherent differences in the 
input physics, most importantly: 

• The initial-final mass relation for white dwarfs; 

• For which systems does mass transfer occur stably, and when is it unstable? 

• When mass transfer is stable, what is the mass transfer rate? 

• How much mass (and angular momentum) is lost in the wind of stars? 

• The evolution of hydrogen-poor helium-burning stars is very important for SNIa pro- 
genitors. Not much is known yet about these objects regarding e.g. the stellar evolution, 
winds or mass transfer stability. 
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